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Summary 


The  research  supported  under  Contract  F49620-92-C-0033  represents  an  effort  to  un¬ 
derstand,  through  observational,  numerical,  and  theoretical  means,  the  role  of  unbalanced 
motions  such  as  gravity  waves,  inertia-gravity  waves,  equatorially  trapped  waves  and  in¬ 
ertial  instabilities  in  the  earth’s  atmosphere,  and  to  develop  suitable  parameterizations  of 
their  effects.  These  waves  and  instabilities  redistribute  momentum  vertically  and/or  hori¬ 
zontally  and  irreversibly  alter  the  distribution  of  momentum,  heat  and  trace  constituents 
when  undergoing  turbulent  breakdown  and  absorption. 

This  report  continues  our  investigation  with  a  theoretical  study  of  equatorial  inertial 
instability.  In  addition  to  the  materia!  presented  herein,  we  are  currently  engaged  in  studies 
of  gravity  and  equatorial  wave  excitation  by  latent  heating  in  the  tropical  atmosphere,  and 
the  role  of  gravity  wave  drag  in  the  bifurcation  properties  of  the  polar  vortex  in  the  winter 
middle  atmosphere.  These  areas  are  being  investigated  in  collaboration  with  Drs.  Francis 
Crum  and  Donal  O’Sullivan  of  NWRA,  who  have  been  partially  supported  by  the  AFOSR 
Contract.  The  substance  of  these  investigations  will  be  reported  later. 
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1.  Introduction 


Inertial  instability  may  arise  in  conservative  axisymmetric  flow  near  the  equator  when 
there  is  nonzero  latitudinal  shear  7  =  u„.  In  this  case  potential  vorticity  (PV)  on  one 
side  of  the  equator  is  anomalous,  satisfying  a  necessary  condition  for  centrifugal  parcel 
instability  (Dunkerton,  1981).  Analytic  eigenmodes  of  symmetric  instability  in  zonally 
uniform  shear  on  an  equatorial  beta-plane  were  obtained  by  Dunkerton  and  independently 
by  Stevens  (1983).  Subsequently  Boyd  and  Christidis  (1982)  and  Dunkerton  (1983)  found 
that  low-wavenumber,  zonally  nonsymmetric  instabilities  have  larger  growth  rates  than 
synunetric  instability  when  (a)  the  vertical  wavenumber  m  is  near  or  below  a  neutral 
point  of  symmetric  instability 

H  -  .1/2  _  M  (ID 

Pf  -  ^neut  -  ^2 

{N  is  static  stability  and  0  =  2Cl/ a  —  2.29-10"^^m“^s~^)  and  (b)  there  exist  nonzero 
integer  zonal  wavenumbers  s  =  ka  below  a  short-wave  cutoff  fcmw  =  0/'1‘  Although 
symmetric  instability  has  maximum  inviscid  growth  approaching  ■7/2  as  |m|  -♦  00,  the 
nonsymmetric  mode  is  preferred  when  second-order  viscosity  and  diffusion  are  added  and 
conditions  (a,b)  are  met  (Dunkerton,  1983).  Similar  conclusions  pertain  to  midlatitude 
flow  containing  nonsymmetric  inertial  and  barotropic  instabilities  as  shown  by  Stevens 
and  Ciesielski  (1986). 


Anomalous  PV  exists  in  the  tropical  winter  mesosphere  (Dunkerton,  1981;  Hitchman 
and  Leovy,  1986),  upper  troposphere  between  South  America  and  Indian  Ocean  (Liebmann, 
1987),  lower  troposphere  near  Indonesia  (Krishnamurti  et  al,  1985,  1988),  and  Pacific 
south  equatorial  current  (Philander,  1989,  p.  62).  Inertial  instabilities  were  found  in  a 
nearly  inviscid  2D  axisymmetric  model  of  the  troposphere  (Held  and  Hou,  1980)  and  in  3D 
middle  atmosphere  GCMs  (Hunt,  1981;  O’Sullivan  and  Hitchman,  1992).  This  instability 
was  suggested  as  a  possible  explanation  of  layered  structures  in  the  tropical  mesosphere 
(Hitchman  et  al,  1987;  Fritts  et  al,  1992)  and  mesoscale  anomalies  in  the  midlatitude 
upper  troposphere  (Ciesielski  et  at,  1989).  Conditional  ‘synunetric*  instability  is  thought 
to  expliun  frontal  runbands  according  to  Bennetts  and  Hoskins  (1979)  and  several  others. 
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Unfortunately  for  the  theory,  none  of  the  regions  of  anomalous  PV  in  atmosphere  or 
ocean  are  zonally  symmetric.  Observations  and  numerical  models  suggest,  for  example, 
that  inertial  instabilities  in  the  mesosphere  coincide  with  tropical  penetration  of  planetary 
Rossby  waves  from  midlatitudes  (Hitchman  et  al,  1987;  O’Sullivzm  and  Hitchman,  1992). 
For  the  theory  to  be  relevant,  it  mxist  be  generalized  to  a  nonparallel  or  zonally  nonuniform 
basic  state. 

The  purpose  of  this  paper  is  to  develop  a  simple  theoretical  formalism  and  thereby  to 
document  two  types  of  inertial  instability  in  nonparallel  flow:  ‘local’  instability  (stationary 
and  zonally  propagating)  and  ‘global’  instability  (symmetric  and  nonsymmetric).  Of  these, 
the  local  instabilities  are  new  and  probably  most  relevant  to  middle  atmosphere  models  in 
which  breaking  planetary  waves  produce  regions  of  signiflcantly  anomalous  PV  near  the 
equator. 

The  theoretical  model  is  described  in  Section  2  including  a  brief  review  of  symmetric 
and  nonsymmetric  instability.  Growth  rates  and  structure  of  nonparallel  instabilities  are 
illustrated  in  Section  3  for  a  wide  range  of  model  parameters  using  a  simple  analytic  basic 
state  resembling  the  onset  of  planetary  wavebreaking.  The  theoretical  interpretation  of 
local  instability  is  discussed  further  in  Section  4  and  the  Appendix. 
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2.  Background 


To  isolate  the  effect  of  horizontal  basic  state  variations,  a  simple  two-dimensional  (x- 
y)  model  is  desirable,  analogous  to  the  shallow-water  system.  In  a  vertically  continuous 
atmosphere,  these  equations  describe  perturbations  oscillating  in  height  with  constant  si¬ 
nusoidal  variation.  Vertical  wavelength  is  then  contained  in  Lamb’s  parameter  e  (Andrews 
et  al,  1987)  and  is  specified  a  priori.  In  this  section,  the  simple  model  is  developed,  its 
solution  method  described,  and  a  brief  review  given  of  parallel  flow  instabilities  (several 
features  of  which  carry  over  to  the  nonparallel  case). 

a.  Theoretical  model 

The  hydrostatic  primitive  equations  on  an  equatorial  beta-plane  are 


tt(  +  ttu*  -1-  u(tty  -  fiy)  -f-  wu,  -b  ^*  =  X 

(2.1a) 

vt  +  tt(v,  +  /?y)  +  uuy  +  tuv,  +  =  y 

(2.1b) 

<f>tt  -f  u^tz  +  u^,y  -b  u;(JV*  -b  <^,,)  =  Q 

(2.1c) 

u*  +  Vy  +  —  (pow).  =  0 

Po 

(2.1d) 

where  u,v,u;  are  zonal,  meridional,  and  vertical  velocity,  ^  is  geopotential,  Po  is  basic  state 
density  p,exp{—zfH),  and  X,  Y,  Q  are  dissipative  terms.  Eqs.  (2.1)  are  essentially  those 
of  Andrews  and  McIntyre  (1976)  written  in  log-pr^sure  notation  of  Holton  (1975). 

To  formulate  an  idealized  model  of  instability  the  dependent  variables  may  be  expanded 
as 

tt(x,y,z,t)  =  tf(i,y)-ftt'(x,y,z,t)  (2.2) 

and  similarly  for  v,w,d>-  U{x,y)  is  a  barotropic  basic  state  independent  of  height  and  time 
for  the  purpose  of  linear  stability  analysis.  Linearized  perturbation  equations  are  then 

uj  -f-  flu',  +  w'(fly  -  /?y)  +  ^',  =  A  +  KV]ju'  -  urn^u' 

v[  -I-  flu'  -1-  /?yu'  +  t^y  =  B  +  KVjjv'  -  t/m*u' 

4/,  +  +  e-*(ui  +  ry  =  C  +  -  KmV' 
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(2.3a) 

(2.3b) 

(2.3c) 


under  the  following  assumptions: 


1)  The  basic  state  velocity  field  is  barotropic,  temporally  constant,  and  horizontal  {W  =  0). 
More  generally  we  could  regzu'd  the  mean  flow  as  slowly  varying  in  height  and  time  (Boyd, 
1978). 

2)  Advection  of  basic-state  temperature  is  neglected,  i.e.  u'  •  =  0.  This  term  has  a 

counterpart  u'  •  normally  retained  in  the  shallow- water  system.^ 

3)  Pertmbations  vary  as  exptmz,  and  m*  »  whereupon  e  «  rn^/N^  as  in  (1.1). 

4)  |0„|  <  N\ 

5)  Dissipative  terms  are  written  as  second-order  viscosity  and  diffusion. 


In  (2.3a-c)  the  zonally  symmetric  c<>mponent  of  U  was  kept  on  the  Ihs,  i.e. 


U  -*xi{y)  +  U{x,y)  (2.4a) 

V-  V{x,y)  (2.4b) 

and  on  the  rhs 

-A  =  (Cfu;  +  Vu'J  +  (tt'Cf,  +  v’Uy)  (2.5a) 

-B  =  {Uv',  +  Vv[)  +  (u'n  +  vX)  (2.5b) 

-C  =  {U<i/,  +  V<f>’^)  (2.5c) 


h.  Solution  method 

To  solve  (2.3a^c)  a  semi-implicit  method  was  used  in  which  variables  were  expanded  in 
zonal  harmonics  (letting  k  be  integers  for  the  moment) 

^In  real  atmoapheres  the  bask  state  will  include  variations  of  static  stability  induced  by  potential  vortkity 
anomalies  (Hoskins  et  al,  1985).  For  simplkity,  only  the  variation  of  induced  horisontal  velocity  U,  V  was 
retained  in  the  definition  oi  A,  B,  C.  This  assumption  b  reasonable  in  the  tropkal  middle  atmosphere  where 
iaeident  Rossby  waves  have  approximately  horisontal  ray  paths. 
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tt 


ikx 
exp - 


(2.6a) 


Ak 

Bi 

Ct 


Terms  on  the  Ihs  of  (2.3a-c)  were  treated  implicitly  in  time,  i.e. 

=  “(«fc 

so  that  for  each  harmonic 


-wilt  +  5k(fi,  -  Py)  +  =  -*AJ  -  wuj 

+(Dt/k  +  ^yHt  + 

ay 


where  u  s  u;  —  kU/a,  u>  =  2i/  At,  and 

A*  =  A*  +  diffusion  terms,  etc. 

After  some  manipulation  (2.8a-c)  reduce  to  a  geopotential  equation 

for  each  harmonic,  where 

A  s  0y{Py  -  7)  “  w*- 


(2.6b) 

(2.6c) 

(2.7) 

(2.8a) 

(2.8b) 

(2.8c) 

(2.9) 

(2.10) 

(2.11) 


D,  whose  definition  is  omitted  for  brevity,  includes  the  explicit  forcing  and  solutions  from 
previous  time  step.  Eq.  (2.10)  was  solved  by  discretizing  in  y  with  64  grid  points  and 
using  a  tridiagonal  algorithm.  The  domain  was  periodic  in  x  and  extended  from  -4y,  to 
+6y,  where  y,  =  '7/2/?.  There  are  no  critical  latitude  (w  =  0)  or  inertial  latitude  (A  = 
0)  singularities  when  Cj  is  complex.  A  spectral  rather  than  pseudospectral  method  was 
used  to  evaluate  D  since  in  this  case  Ut,  Vi  =  0  for  any  |fc|  1,  that  is,  a  purely  wave  1 
contribution  to  basic  state  variation  in  z  was  assumed.  Product  terms  At,  etc.  could  be 
determined  efficiently  without  transforms.  In  the  numerical  integrations 

«  =  1f(y  -  yc)  +  fio  (2.12) 

where  q  auid  uq  are  constants.  Fields  of  t/,  V  are  specified  in  Section  3. 


c.  Instability  of  uniform  parallel  flow 

When  U  —  V  =0,  harmonic  perturbations  a  exp  i{kx  —  wt)  satisfy 


-Cju  +  v(7  -  0y)  +k<t>  =  0 

(2.13a) 

+WU  +  Pyu  +<f>y  =0 

(2.13b) 

—Cjt<f>  +  A;tt  +  Vy  =  0 

(2.13c) 

where  u>  =  uj  —  ku.  The  geopotential  equation  derived  from  (2.13a-c)  is  just  (2.10)  with 
P  =  0  and  u>  redefined  as  intrinsic  frequency.  This  equation  is  linear  in  the  ‘eigenvalue’  e. 
Boyd  (1978)  derived  the  meridional  velocity  equation 


2£wfc'7 
cw’  —  k^ 


=  0 


(2.14) 


which  is  nonlinear  in  £  unless,  for  example,  >  A;’  in  which  case 


2k'i 

«vv  +  — -  « 


-(/?-ir»)  +  cA 
w 


«0 


lfib  =  0 


v„  -  tv  [Py[Py  -  7)  -  w*]  =  0 


(2.15) 


(2.16) 
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and  the  eigencondition  for  symmetric  instability  is 


w*  +  t/4  = 


(2n  +  1)NI3 


(2.17) 


(Dunkerton,  1981).  Solutions  for  v  are  Hermite  functions  (polynomials  times  a  Gaussism) 
centered  about  a  ‘shifted  equator’  y,  =  'y/2)3,  the  center  of  zmomalous  vorticity  region. 
For  n  =  0  the  neutral  point  of  symmetric  instability  (hereafter,  ‘neutral  point’)  is  given  by 
(1.1).  For  small  k  and  =  0, 


,  NI3  k0 

m  W6 


At  6  —  Cn«uti 


(jJ^  « 


(2.18a) 


(2.18b) 


Eqs.  (2.18a,b)  were  derived  in  nondimensional  form  by  Boyd  and  Christidis  (1982). 

At  finite  k  the  eigenproblem  (2.10)  or  (2.14)  must  be  solved  numerically.  Growth  rates 
for  'y  =  10“®s~^  are  shown  in  Fig.  1  illustrating  the  relative  importance  of  symmetric 
and  nonsymmetric  modes.  Curves  were  obtained  from  the  geopotential  equation  by  a 
shooting  method;  symbols  correspond  to  values  from  the  time-dependent  model.  They  are 
in  excellent  agreement  except  at  very  small  growth  rates,  where  the  time-dependent  model 
could  not  determine  growth  rate  accurately,  and  at  large  e,  where  latitudinal  resolution 
was  inadequate. 

An  approximate  model  derived  from  (2.15)  contains  eigenfrequencies  similar  to  the 
exact  values  of  Fig.  1  over  a  large  portion  of  parameter  space.  This  model  begins  with 
constant  qr  and  assumes  that  y-variations  of  intrinsic  frequency  can  be  neglected  outside 


■v„  —  V 


+  €A  «  0 


(2.19) 


where  a;o  is  a  constant  (complex)  intrinsic  frequency  within  the  unstable  region.  It  follows 


“.r  +  V  =  MJi  -  +  M 

1  -  4^  lm|  ^  '  WqC 


(2.20) 
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GROWTH  RATE  W;  /y 


ZONAL  WAVENUMBER  k 


Fig.  1:  Growth  rate  as  a  fonction  of  zonal  wavenumber  in  parallel  flow,  for  7  s= 

e  in  units  m~’s*.  Symbols  denote  results  of  time-dependent  numerical  model;  solid 
and  dashed  curves  show  analytic  results  from  a  shooting  method.  (The  dashed  curve 
had  no  symmetric  instability  because  e  <  Cmu*-) 
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where 


Wo  =  w  —  k'^ 


2  /3  —  2wfc 


(2.21) 


- -  '  2{;9>  -  kW)  '  ' 

A  formal  jiistification  of  approximations  leading  to  (2.20)  will  not  be  given  except  to  note 
that  terms  of  arise  from  A  at  large  y  and  cannot  be  ignored  for  latitudinal  trapping 

and  quantization  of  eigenfunctions  near  the  short-wave  cutoff.  The  approximate  dispersion 
relation  (2.20)  has  the  desired  behavior  at  small  and  large  k.  In  between,  agreement  with 
exact  growth  rates  and  phase  speeds  is  good  for  £  >  Cncut  (not  shown).  At  smaller  £  the 
cubic  equation  fails  in  two  ways:  growth  rates  are  overestimated  at  intermediate  k  and 
there  is  no  Kelvin  wave  instability  at  small  k  below  the  neutral  point  (Boyd  and  Christidis, 


1982). 
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3.  Instability  of  nonparallel  flow 

The  time-dependent  model  of  Section  2  was  used  to  determine  the  growth  of  unstable 
modes  in  nonpaiallel  flow.  Some  characteristics  of  these  modes  are  now  described.  The 
basic  state,  independent  of  height  and  time,  was  designed  as  a  simple  analog  of  incipient 
planetary  wavebreaking  in  the  tropical  mesosphere.  The  evolution  of  inertial  instability 
in  a  time-dependent  Rossby-wave  critical  layer  was  recently  discussed  by  O’Sullivan  and 
Hitchman  (1992)  in  a  three-dimensional  model.  In  order  to  understamd  the  selection  of 
unstable  modes  it  will  be  worthwhile  to  explore  the  parameter  dependence  of  inertial  in¬ 
stability  in  a  simpler  two-dimensional  model. 


a.  Basic  state 

The  basic  state  was  defined  by 

1  +  tanh  cos(x/o) -h  (3.1) 

In  (3.1),  the  latitude  of  maximum  shear  is  ye,  ^  is  dimensionless  amplitude  of  basic  state 
variation,  y,  =  ir/2)9  is  the  center  of  anomalous  zonal-mean  vorticity  as  defined  previously, 
and  So  is  a  constant  mean  flow  independent  of  x,  y.  The  profile  (3.1)  is  shoMm  in  Fig.  2  for 
tfe  =  1.8y,,  t/f  =  1,  and  So  =  0. 

In  Fig.  2,  latitudinal  shear  at  y  =  ye  was  increased  (decreased)  relative  to  the  zonal 
mean  by  50%  at  180  (0)  degrees  longitude.  The  flow  was  inertizdly  unstable  at  all  x,  but 
more  so  at  the  center  of  figure.  This  region  could  represent  a  zone  of  enhanced  latitudinal 
shear  set  up  by  Rossby  waves  penetrating  from  latitudes  north  of  the  equator.  In  reality 
there  would  be  some  latitudinal  phase  tilt  during  critical  layer  development  (O’Sullivan 
and  Hitchman,  1992)  -  an  unnecessary  complication  for  stability  analysis. 

The  meridional  component  of  basic  state  V(x,y)  was  included  in  all  simulations,  i.e., 

Cf,  +  V;  =  0.  (3.2) 


fi(y)  +  u{x,  y)  =  7(y  -  ye)  - 
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LONGITUDE 

2:  Profile  of  nonparallel  basic  state  (3.1)  with  0  =  1,  Contours  show  dimensionless 

parcel  displacement  growth  rate  iZ/^.  Contour  interval  0.1,  with  0.1  and  0.4  contours 
omitted. 


This  had  a  minor  effect  compzired  to  the  zonal  component  (3.1).  Although  the  details  are 
omitted,  it  can  be  stated  that  all  of  the  results  shown  here  were  similar  when  V  =  0.  Local 
instabilities  were  also  found  when  U  =  0  and  V  =  ^(x)  (localization  due  to  dVjdx).  The 
effect  of  V  was  generally  to  disrupt  rather  than  enhance  the  instability  in  comparison  to 
parallel  flow.  These  cases  were  felt  to  be  less  interesting  and  were  not  investigated  further. 


6.  Latitude  of  maximum  shear 

The  latitude  of  maximum  latitudinal  shear  was  j/e  according  to  (3.1).  On  the  other 
hand,  the  maximum  rate  of  ‘parcel  instability’  \s  R  =  +  ?)  where  f  is  relative 

vorticity.  This  quantity  is  contoured  in  Fig.  2;  it  did  not  always  maximize  at  y  =  ye  but 
was  a  function  of  ye  (and  other  parzuneters),  having  (for  0  =  1)  an  overall  maximum  value 
R  —  0.75^  when  y  =  =  l.Sy*.  Maximum  R  was  located  north  of  y,  because  the  wave’s 

contribution  to  basic  state  expanded  the  zone  of  instability  at  180®  (recall  that  y,  was  the 
center  of  anomalous  zonal-mean  vorticity).  However,  as  |yel  — ^  oo  the  relative  vorticity  at 
|ye(  became  small  relative  to  jSy,  so  the  region  of  anomalous  vorticity  next  to  the  equator 
returned  to  a  parallel  configuration.  For  parallel  flow  the  maximum  =  7/2  at  y  =  y,. 

Growth  rates  of  the  most  unstable  mode  as  a  function  of  ye  are  shown  in  Fig.  3  for 
7  =  4*  10~®  s“',  e~^  =  256  m*s~*,  0  =  1,  and  uq  =  0.  For  these  calculations  8  zonal 
harmonics  were  used,  and  At  =  900  s.  Vertical  wavelength  ~5  km  for  this  choice  of  e 
when  N  =  .02s"  The  asymptotic  form  of  mode  emd  growth  rate  were  realized  after 
a  few  days  in  most  cases,  although  as  will  be  seen  later,  this  time  scale  depended  on 

7.  The  results  appear  reasonable:  growth  rate  peaked  at  y*  «  1.7y,  and  for  large  (yd 
instabilities  developed  as  in  parallel  flow  (not  shown).  However,  the  interpretation  of 
Fig.  3  is  slightly  more  complicated,  because:  (l)  Instabilities  near  the  center  of  the  plot 
were  exactly  stationary  (‘s’)  and  grew  about  twice  as  fast  as  in  parallel  flow  -  a  larger 
variation  than  expected  from  R  alone.  Peak  growth  occurred  a  little  to  the  right  of  1.5y,. 
(2)  Instabilities  adjacent  to  those  labeled  ‘s’  were  zonally  propagating  but  not  global  in 
extent. 
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The  basic  state  zonal  flow  [u(j/)  +  U{x,y)  in  (3.1)]  at  location  of  maximum  ^(180®,!/) 
also  varied  with  ye',  it  was  zero  at  this  point  when  ye  «  2.2y«.  This  variation  affected  the 
growth  rate  in  addition  to  the  variation  of  R.  Growth  was  optimized  for  a  combination 
that  maximized  R  and  minimized  the  m  situ  mean  flow.  The  dependence  of  growth  rate  on 
mean  flow  speed,  or  breakdown  of  Galilean  invariance,  is  characteristic  of  ‘local’  instability 
of  nonparallel  flow  as  disciissed,  e.g.,  by  Pierrehumbert  (1984).  We  return  to  this  point  in 
the  next  subsection  and  in  Section  4. 

Instability  structure  for  =  1.8y,  is  shown  in  Fig.  4.  The  disturbance  was  exactly 
stationary  and  confined  to  the  middle  half  of  the  domain.  Velocity  was  directed  from  low 
to  high  pressure  indicating  a  dynamical  source  of  instability  (Coriolis  force).  Curiously 
the  geopotential  phase  tilted  slightly  although  no  tilt  was  imposed  on  the  basic  state  (3.1). 
Identical  results  were  obtained  with  doubled  resolution  in  x,  y  and  At  =  450  s. 

Fig.  4  exemplifies  ‘local  stationary’  instability.  Note,  (1)  although  localized  in  x  it  is 
not  infinitesimal  in  size.  This  indicates  that  terms  in  addition  to  v'Uy  are  important  in  the 
dynamics,  as  expected  from  Section  2c.  An  oversimplified  model  could  be  formulated  in 
which  u'C/y  was  the  only  term  in  A,  all  linear  terms  oc  dfdx  were  dropped,  and  B,C  =  0. 
Growth  rate  in  this  case  could  be  obtained  from  the  symmetric  formula  (2.18a  with  fc  =  0) 
but  was  a  function  of  x,  implying  a  5-function  catastrophe.  (This  behavior  was  verified 
numerically.)  (2)  Absence  of  zonal  propagation  over  a  finite  rzmge  of  parameters  (in  this 
case,  ye)  apparently  requires  the  zonally  symmetric  (s  =  0)  component  of  perturbation, 
an  important  part  of  the  total  wave  field.  Most  of  the  perturbation  energy  was  contained 
in  the  lowest  four  harmonics.  When  the  s  =  0  component  was  artificially  excluded  from 
the  numerical  model,  real  phase  speed  was  nonzero  except  where  its  trajectory  crossed  the 
Cr  =  0  axis  at  a  point  near  ye  =  1.8y,. 

In  parallel  flow,  nonsymmetric  instability  propagates  zonally  except  in  the  special  case 
er{k,€)  =  —Ho  (Dunkerton,  1983;  Stevens  and  Ciesielski,  1986).  Only  symmetric  instabil¬ 
ity  is  trivially  ‘stationary’  regardless  of  parameters.  In  nonparallel  flow,  local  stationary 
instability  is  possible.  It  is  in  a  sense  ‘locally  symmetric’  although  this  description  ignores 
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4:  Structure  of  local  stationary  mode  for  -7  =  4  •  =  256m*s“ 

flo  =  0,  and  j/a  =  l.Sya.  Negative  ^-contours  are  dashed. 


the  spatial  variation  of  instability  structure  (e.g.,  Fig.  4). 


c.  Dtpendtnee  on  So 

The  term  Sq  in  (3.1)  added  a  constant  mean  flow  independent  of  x,y.  In  parsdlel  flow 
So  would  have  no  effect  other  than  to  shift  phase  speed  by  this  constant.  For  local  modes 
in  nonparallel  flow,  growth  rate  is  altered. 

Fig.  5  shows  a;,-/7  for  the  same  parameters  as  Fig.  2  but  with  Pe  =  1.8y,  and  So  varied 
over  a  large  range.  Local  stationary  instability  occurred  near  the  point  of  zero  mean  flow, 
but  for  large  |So|,  the  instability  was  displaced  off  center  and  forced  to  propagate  in  the 
direction  of  So. 

An  example  of  this  zonally  propagating  or  ‘advected’  instability  is  shown  in  Fig.  6  for 
So  =  -30  ms~^.  Propagation  was  to  the  west,  and  the  pattern  was  observed  to  repeat 
every  six  days.  Obviously  this  was  insufflcient  time  for  global  traverse,  and  in  any  case, 
the  disturbance  re-entering  &om  the  right  had  little  role  to  play  in  the  sequence. 

It  is  generally  thought  that  local  instabilities  belong  to  a  class  defined  by  ‘absolute’ 
instability  (Pierrehumbert,  1984).  Absolute  instability  is  sensitive  to  a  mean  flow  that 
displaces  the  center  of  mode  away  from  the  source  of  instability  and  thereby  lowers  the 
effective  growth  rate.  When  the  mean  flow  is  stable  at  large  |x|,  a  sufficiently  large  |iio| 
can  stabilize  the  problem.  This  did  not  happen  at  the  edges  of  Fig.  5  because  of  the 
zonally  symmetric  component  of  sheau*  and  periodic  boundary  conditions;  advected  modes 
simply  wrappe4  around  and  continued  to  exist  as  ‘global’  modes  in  z,  with  smaller  growth 
rate.  A  sponge  region  near  z  =  0,  to  simulate  an  infinite  domain,  dramatically  altered  the 
advected  mode  as  discussed  at  the  end  of  this  section.  The  interpretation  of  local  modes 
as  absolute  instability  is  deferred  to  the  Appendix. 

Of  the  two  sub- types  of  local  instability  (stationary  and  advected),  the  stationary  va¬ 
riety  seems  more  relevant  to  equatorial  Rossby-wave  critical  layers,  for  instance,  when  a 
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quasi-stationary  Rossby  wave  induces  a  region  of  significantly  auiomalous  potential  vortic- 
ity  near  the  zero-wind  line.  Here  the  advecting  mean  flow  is  small.  That  is  presumably 
why  instabilities  simulated  by  O’Sullivan  and  Hitchman  (1992)  were  quasi-stationary  (‘lo¬ 
cally  symmetric’)  and  followed  the  anomalous  PV  during  critical  layer  development. 


d.  Transition  to  longitudinally  ‘global’  instability 

The  value  of  tj;  in  (3.1)  determined  the  zonal  variation  of  shear.  It  was  found  to  no 
surprise  that  as  ^  decreased  there  was  a  transition  to  global  instability,  as  shown  in  Fig.  7 
(‘p’  signifying  a  globally  propagating  or  parallel-flow  instability).  Growth  rate  was  reduced 
as  0  — 0.  Parameters  were  the  same  as  before,  but  with  ye  =  1.8j/„  So  =  0,  and  variable 
The  structure  became  essentially  that  of  nonsymmetric  instability  in  parallel  flow  (Stevens 
and  Ciesielski,  1986),  with  enhanced  amplitude  near  the  most  unstable  region  (see  Fig.  8, 
for  0  =  0.2).  These  were  global  modes  that  recycled  through  the  entire  domain.  (Following 
Pierrehumbert,  ‘global’  refers  only  to  the  longitudinal  extent  of  eigenmodes.) 

Results  at  very  small  ^  depended  on  initial  conditions.  A  wavenumber  1  initial  distur- 
bsmce,  used  in  most  simulations  here,  produced  wave  1  nonsymmetric  instability.  Initial 
wave  2  led  to  wave  2  instability  at  small  with  slightly  larger  growth  rate.  Wave  3  was 
near  the  short-wave  cutoff  and  less  unstable. 

The  time-dependent  model  simulated  only  the  most  unstable  mode;  however,  we  expect 
local  and  global  modes  of  instability  to  coexist,  at  least  in  cases  such  as  Fig.  2  where  the 
flow  is  unstable  at  all  x.  Whether  global  modes  of  inertial  instability  are  important  in  the 
middle  atmosphere  or  anywhere  else  is  xmcertain.^ 

^There  k  an  intriguing  visible-light  photograph  of  Jupiter  taken  by  the  Hubble  Space  Telescope  on  May  28, 
1991  revealing  horisontaUy-tflted,  banded  structures  next  to  the  equator  (see,  e.g.,  the  cover  page  of  AGU 
publication  Earth  in  Space,  January  1992  issue).  This  would  suggest  horisontally-divergent  circulations 
perhaps  due  to  a  dynamical  instability. 
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Structure  of  global  instability  for  parameters 


e.  Variation  of  c 


The  value  of  in  Figs.  3-8  (256  m*s~*)  was  near  the  neutral  point  =  305  m^s"* 
for  '7  =  4- 10~®  s~^.  Vertical  wavelength  was  ~5  km  when  N  =  .02s~^.  [The  reader  may 
refer  to  (1.1)  for  other  N,  and  recall  the  transformation  of  variables  in  Dunkerton  (1983) 
allowing  generalization  to  arbitrary  values  of  shear.] 

Recalling  the  discussion  of  Section  2c,  the  ‘neutral  point’  limits  symmetric  instability 
in  parallel  flow.  There  is  no  neutral  curve  at  small  e  (except  at  small  k  adjacent  to  the 
neutral  point;  see  Boyd  and  Christidis,  1982)  to  prevent  instability  at  nonzero  s  =  ka.  For 
7  =  4- 10“®  s“^,  several  unstable  integer  wavenumbers  are  allowed  below  the  short-wave 
cutoff  (s  =  0-3).  Growth  rate  varies  monotonically,  increasing  with  £  to  a  limiting  value 
7/2  (cf.  Fig.  1). 

In  nonparallel  flow,  the  growth  of  local  instability  depended  on  £  in  a  similar  way  as 
shown  in  Fig.  9.  The  same  parameters  were  used  in  (3.1)  with  =  1.8y,,  =  1,  and 

Qo  =  0.  Growth  rate  varied  monotonically,  reaching  a  plateau  at  higher  £.  Over  most  of  this 
range  the  instability  was  stationary,  displaying  slow  eastward  propagation  only  at  small  £. 
It  seems  clear  from  this  example  that  nonparallel  instability  exists  below  the  neutral  point 
whether  one  adopts  a  ‘zonal  mean-flow’  value  (f^ut  ~  305  m*s“*)  or  hypothetical  ‘local’ 
value  (fnJIit  ~  1545  m*s”*)  at  center  of  domain. 

The  structure  of  instability  at  small  £  had  several  interesting  features  (not  shown):  the 
latitudinal  scale  was  enlarged,  as  in  classical  equatorial  wave  theory,  so  instability  extended 
into  the  stable  part  of  the  domain.  Its  structure  resembled  a  Rossby  wave  in  the  northern 
hemisphere  and  eastward  inertia-gravity  wave  in  the  southern  hemisphere  (Matsimo,  1966). 
The  real  phase  speed,  though  slow  eastward,  was  for  practical  purposes  stationary  with 
respect  to  the  ground.  The  disturbance  propagated,  with  respect  to  the  fluid,  westward 
north  of  the  equator  and  eastward  south  of  the  equator.  Because  of  latitudinal  shear,  a 
fused  Rossby /inertia-gravity  structure  could  be  excited  by  equatorial  instability  that  acted 
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as  a  kind  of  ‘wavemaker’  in  this  instance.^  The  instability  was,  however,  local  in  x. 

The  structure  of  instability  at  large  e  resembled  Fig.  4  but  was  contracted  in  latitude, 
remaining  local  in  x. 

Fig.  9  suggests  that  as  in  parallel  flow,  scale  selection  at  finite  |m|  is  unrealized  in 
the  simple  model.  Dunkerton  (1981,  1983)  argued  that  vertical  diffusion  selects  finite 
|m|.  With  diffusion  there  is  a  net  growth  rate  =  w,-  —  utN^.  In  the  real  world  it  is 
uncertain  whether  diffusion  causes  scale  selection  or  is  a  posteriori  ca\ised  by  instability. 
(In  the  mesosphere,  breaking  internal  gravity  waves  will  help  stabilize  the  flow.)  Other 
factors  may  be  important:  (1)  Stability  depends  on  basic  state  transience;  instabilities 
must  grow  rapidly  compzu^ed  to  mean-flow  change.  In  the  tropical  mesosphere,  planetary 
waves  will  modify  the  basic  state  on  a  time  scale  of  1-5  days.  (A  Lagrangian  time  scale  is 
what  matters  inasmuch  as  the  basic  state  conserves  potential  vorticity  on  fluid  parcels  and 
inertial  instability  is  to  some  extent  a  ‘parcel’  instability  dependent  on  anomalous  PV.) 
It  would  therefore  be  consistent  to  require  7  >  1-2*10"®  s~^  for  instability  (cf.  Hitchman 
and  Leovy,  1986;  Hitchman  et  al,  1987).  (2)  Realistic  basic  states  contain  a  continuous 
spectrum  of  horizontal  and  vertical  motions.  The  initial  evolution  need  not  be  dominated 
by  the  most  unstable  mode,  in  fact  the  instability  at  finite  amplitude  need  not  be  modal. 
(3)  The  meridional  circulation  itself  has  a  weak  stabilizing  influence  at  large  e  (Section  3f). 

/.  Effect  of  meridional  circulation 

Cross-equatorial  shear  in  the  middle  atmosphere  is  due  to  diabatic  advection  by  a  mean 
meridional  circulation.  Including  in  (2.3)  a  zonally  symmetric  component  9o>  independent 
of  latitude,  had  a  minor  effect  on  local  stationary  instability  as  shown  in  Fig.  10.  For  this 
series  7  =  4*  10”*  s“*,  c”*  =  256  m*s"*,  y*  =  1.8y„  ^  =  1,  and  Cq  =  0-  The  location  of 
instability  was  shifted  north  and  east  when  i^o  >  0  (not  shown).  Growth  rate  was  slightly 
faster  in  most  cases  with  €~^  >  256  m*s~*  plotted  in  Figs.  3,5, 7,9  when  Oq  =  5  ms”^.  At 

’A’ontal  excitation  of  inertia-gravity  waves  in  a  localised  region  of  conditional  ‘symmetric’  instability  was 
noted  by  Thorpe  and  Rotnnno  (1989)  and  Jones  and  Thorpe  (1992). 
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large  e  there  was  a  reduction  of  growth  suggesting  the  possibility  of  weak  scale  selection 
when  So  >  0.  The  physical  explanation  may  be  that  So  advects  the  instability  away  from 
maximum  R,  the  effect  being  more  dramatic  for  high-£  modes  with  short  meridional  scale. 


g.  Results  for  weak  shear 

The  preceding  results  raise  two  questions  about  instability  in  nonparallel  flow,  (l)  Ob¬ 
serve  from  Figs.  4,6  that  local  instability  had  the  same  zonal  scale  as  basic  state  variation. 
There  was  at  most  one  stationary  or  two  propagating  ceils  in  the  x-direction.  In  general, 
we  expect  multi-cell  nonsymmetric  instability  in  ‘locally  parallel*  flow.  Under  what  circum¬ 
stances,  then,  can  a  wavetrain  of  many  ceils  develop  in  the  unstable  part  of  the  flow?  (2) 
Advected  instability  was  observed  at  large  |uo|  because  the  flow  was  unstable  at  all  x  and 
periodic  boundary  conditions  were  assumed.  Real  flows  could  be  stable  at  some  longitude, 
fa  advected  instability  possible  in  this  case?  The  first  question  will  be  addressed  now  and 
the  second  in  the  next  subsection. 

It  was  seen  in  parallel  flow  that  the  minimum  zonal  scale  for  latitudinal  trapping 
of  eigenfvinctions  (decay  at  large  |y()  is  defined  by  a  short-wave  cutoff  kmaxl  =  The 
argument  could  be  extended  to  nonparallel  flow;  minimum  zonal  scale  depends  on  shear.  In 
the  case  illustrated  above,  it  was  impossible  to  accomodate  several  cells  in  the  x-direction. 
For  multiple  cells  7  must  be  reduced.^ 

Fig.  11  shows  an  example  in  which  7  =  0.5-10~®  s~^,  c  =  4  m~  V,  Pc  =  1.8y„  ^  =  1,  and 
Go  =  0.  Numerical  integration  \ised  32  harmonics  (note  change  of  axes) ;  results  are  shown 
at  240  days.  For  this  value  of  7  a  long  time  was  required  to  obtain  the  mode,  due  to  the 
assumed  7  and  wave  1  initial  condition.  (Less  time  would  have  been  required  from  high-k 
initial  conditions.)  This  example  demonstrated  that  multi-cell  instability  is  possible  when 
the  characteristic  scale  of  mean-flow  variation  is  much  greater  than  k~l,  =  7/^.  Another 

^That  is,  in  the  configuration  (3.1).  There  may  exist  situations  where  7  is  unimportant  outside  the  region 
of  aaomalons  PV,  such  that  trapping  of  eigenfunctions  is  ensured  by  the  basic  state  alone,  ^art  fi-om  k. 
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Structure  of  multi-cell  instability  in  weak  shear  7  =  0.5  •  10~ 


example  of  small  7  is  discussed  in  the  Appendix. 

What  about  the  opposite  situation,  when  the  unstable  region  is  smaller  than  A;~L? 
Profiles  of  the  form  (3.1)  could  be  generalized  in  order  to  confine  the  most  unstable  flow: 

—  cos(i/a)  =  2sin*(i/2a)  —  1  — v  2sin'‘(i/2o)  —  1  (3.3) 

Setting  n  =  16,  for  example,  the  region  of  most  unstable  flow  was  less  than  60“  wide. 
Results  indicated  that,  for  parameters  of  Fig.  4,  the  local  stationary  mode  was  essentially 
unchanged  by  reducing  the  size  of  unstable  region,  although  growth  rate  was  reduced 
slightly  (by  about  25%  when  n  =  16).  (An  explicit  version  of  the  time-dependent  code  was 
Tised  to  obtain  this  result  and  to  validate  other  results  of  the  semi-implicit  model.)  This 
supports  its  interpretation  as  a  local  mode  that  depends,  not  on  global  average  stability, 
but  local  stability  within  the  most  unstable  part  of  the  domain  (Pierrehumbert,  1984) . 

To  be  stire,  the  concept  of  ‘local  stability’  is  imprecise  when  (as  in  this  case)  instability 
and  basic  state  share  the  same  zonal  scale.  WKB  analysis  like  that  of  Pierrehumbert  (1984) 
is  formally  valid  only  in  slowly  varying  mean  flow.  Perhaps  a  more  accurate  measure  of 
stability  would  be  obtained  by  averaging,  say,  over  a  half-wavelength  of  perturbation. 


h.  Sponge  region  in  x 

As  mentioned  in  Section  2c  the  advected  instability  was  altered  by  a  sponge  region 
designed  to  completely  absorb  disturbances  crossing  the  domain  boundary  near  x  =  0. 
Though  artificial,  the  sponge  may  simulate  regions  of  stable  flow  (as  undoubtedly  exist  in 
the  atmosphere)  or  the  tendency  of  instabilities  to  ‘break’  and  dissipate  through  nonlinear 
saturation  before  making  a  complete  circuit  of  the  globe. 

Fig.  12  shows  growth  rate  as  in  Fig.  5  but  with  sponge  region  included.  Instabilities 
were  unaffected  at  small  |&o|  but  were  stabilized  at  large  |6o|-  The  physical  interpretation 
is  that  So  displaced  the  center  of  mode  downstream  from  the  source  of  instability  near  x 
at  180*.  For  sufficiently  large  |ilo|  the  mode  was  pushed  into  the  sponge  region  and  could 
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not  grow.  Evidently  ‘global’  modes  were  stabilized  by  the  sponge  region  also. 
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4.  Discussion 


Results  of  Section  3  documented  local  as  well  as  global  modes  of  inertial  instability  in 
nonparallel  flow.  Local  modes  dominate  inhomogeneous  flow  and  are  likely  more  impor¬ 
tant  in  observed  bzisic  states.  Local  stationary  instabilities  had  the  largest  growth  rate, 
depending  primarily  on  the  local  rather  than  global  stability  of  the  basic  state;  they  oc¬ 
curred  when  (among  other  things)  mean  flow  advection  was  small  in  the  unstable  region. 
L2irge  advection  caused  zonal  phase  propagation,  here  referred  to  as  advected  instability. 
Like  stationary  modes,  advected  modes  had  zero  group  velocity  but  the  center  of  wave 
packet  was  displaced  downstream  from  the  most  unstable  longitude,  and  the  growth  rate 
reduced.  If  uq  was  sufliciently  large  and  the  flow  stable  at  large  |r|,  advected  modes  were 
stabilized. 


Even  for  a  simple  profile  like  (3.1)  the  parameter  dependence  of  inertial  instability  is 
complicated;  examples  of  Section  3  highlighted  behavior  along  a  few  trajectories  in  parame¬ 
ter  space.  One  of  the  more  interesting  results  was  the  dependence  of  local  instability  on  fio, 
illustrating  the  breakdown  of  Galilean  invariance  in  nonparallel  flow.  Zonal  inhomogeneity 
breaks  the  symmetry  of  uniform  parallel  flow  and  introduces  a  new  constant:  the  zonal 
propagation  speed  of  basic-state  pattern.  This  was  zero  in  (3.1);  consequently,  instabilities 
with  zero  group  velocity  remained  in  the  same  location  relative  to  the  mean  flow  pattern 
at  all  times.  (Not  surprisingly  their  growth  maximized  when  the  mode  best  overlapped 
the  most  unstable  region.)  In  WKB  theory,  a  flow  is  said  to  be  absolutely  unstable  when 
such  instabilities  exist,  i.e. 

w<(^)  >  0  (4.1) 


where 


(4.2) 


and  unstable  branches  originate  on  opposite  sides  of  the  real-k  axis  (e.g.,  Pierrehumbert, 
1984, 1986).  The  importance  of  absolute  instability  derives  from  symmetry-breaking,  e.g., 
in  nonparallel/nonuniform  flow,  or  when  there  is  a  ‘wavemaker’  at  some  x  (Huerre  and 
Monkewitz,  1990). 
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Pierrehumbert  (1984)  analyzed  the  two-layer  quasi-geostrophic  model  for  baroclinic  in¬ 
stability  of  nonparallel  flow,  including  a  discussion  of  ‘absolute’  and  ‘convective’  instability 
in  terms  of  local  and  global  modes,  supplemented  by  WKB  analysis.  It  is  unnecessary  to 
duplicate  this  work,  but  simpler  to  quote  three  of  his  main  conclxisions  pertaining  to  locsJ 
modes: 


1)  Growth  rate  of  a  local  mode  is  determined  by  the  maximum  baroelinicity  in  the  domain, 
and  not  by  the  average  baroelinicity. 

This  was  suggested  for  inertially  unstable  flow  (substituting  ‘inertial  instability’  for 
‘baroelinicity’)  by  results  of  Figs.  5,12  demonstrating  the  insensitivity  of  local  stationary 
instability  to  a  sponge  region  near  z  —  0,  and  by  discussion  at  the  end  of  Section  3g  in 
which  growth  rate  was  only  weakly  dependent  on  the  size  of  the  unstable  region.  As  in 
Pierrehumbert  (1984)  the  ‘maximum  baroelinicity’  or  ‘maximum  inertial  instability’  of  the 
basic  state  must  take  into  account  factors  (such  as  boundary  conditions,  pressure  gradients, 
and  So)  in  addition  to  the  parcel  growth  rate  (e.g.,  R).  On  this  account  we  are  reminded 
that  equatorial  inertial  instability  is  not  a  pure  ‘parcel’  instability  (Dunkerton,  1983). 


2)  Vanishing  absolute  growth  rate  at  infinity  is  not  necessary  for  localization. 

This  was  apparent  from  examples  based  on  the  profile  (3.1)  which  was  unstable  at  all 
X  but  displayed  a  preference  for  local  instabilities  in  most  circuiostances. 


3)  The  contrast  between  maximum  and  minimum  baroelinicity  in  the  flow  determines  the 
extent  of  localization  of  the  eigenmodes,  with  high  contrast  favoring  localization. 

bi  a  similar  way  (although  it  was  not  discussed  extensively  above),  the  zonal  scale  of 
local  instability  increased  as  0  was  decreased  below  1,  before  the  onset  of  global  instability. 
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Absolute  instability  does  not  depend  on  the  validity  of  WKB  in  (4.1,2)  since,  by  def¬ 
inition,  the  instability  must  grow  indefinitely  at  a  fixed  point  in  x.  This  property  was 
demonstrated  numerically,  but  not  analytically,  for  local  stationary  instability.  WKB 
analysis  of  inertial  instability  on  an  equatorial  beta-plane  is  difficult  although  suggestive 
of  absolute  instability  as  discussed  in  the  Appendix.  Basic  states  like  Fig.  2  disallow  a  clean 
separation  of  wave  and  mean-flow  scales,  and  a  numerical  demonstration  of  mstability  is 
preferable  in  this  case  (cf.  Jones  and  Thorpe,  1992). 
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5.  Conclusion 


Using  a  simple  model  it  was  shown  that  inertial  instabilities  in  zonally  nonuniform 
cross-equatorial  shear  take  the  form  of  local  and  global  modes.  Local  modes  may  be 
exaKitly  stationary  or  display  zonal  phase  propagation  depending  on,  among  other  things, 
advection  by  the  basic  state.  Both  sub-types  of  local  instability  are  distinguished  from 
global  instability  by  their  zero  group  velocity  and  concentration  of  amplitude  within,  or 
downstream  from,  the  region  of  most  unstable  flow.  These  properties,  including  breakdown 
of  Galilean  invariance,  are  reminiscent  of  ‘absolute’  instability  -  although,  apart  from  the 
Appendix,  none  of  our  results  depended  on  a  WKB  approximation.  Exact  stationarity  of 
phase  over  a  finite  range  of  parameters,  and  measurable  dependence  of  growth  rate  on  the 
size  of  unstable  region,  are  possibly  attributable  to  a  non- WKB  effect. 

It  was  suggested  that  for  equatorial  Rossby-wave  critical  layers  in  the  mesosphere,  local 
stationary  modes  are  important  given  the  strength  of  background  shear  7  and  weakness 
of  m  situ  advecting  current  near  the  zero-wind  line.  Such  a  simplification  is  desirable 
in  this  context  to  avoid  the  complexity  of  global  and  advected  modes.  Existence  of  local 
stationary  modes  locked  to  PV  anomalies  is  consistent  with  middle  atmosphere  simulations 
of  inertial  instability  smd  Rossby  wave  breaking  by  O’Sullivan  and  Hitchman  (1992).  As  a 
follow-up  it  will  be  worthwhile  to  show  the  evolution  of  unstable  modes  in  an  authentic 
Rossby  wave  critical  layer,  and  their  effect  on  potential  vorticity  evolution.  This  will  be 
done  in  the  sequel. 

Local  and  global  modes  of  inertial  instability  share  many  of  the  same  properties,  e.g., 
dependence  of  growth  rate  on  vertical  wavelength,  and  downgradient  momentum  flux. 
Vertical  scale  selection  may  be  achieved  with  diffusion  as  argued  previously,  or  as  the  new 
analysb  suggests,  through  mean  meridional  advection.  Non-modal  instability  (Farrell, 
1982)  on  a  transient  basic  state  will  also  allow  scale  selection  at  finite  amplitude.  It  should 
be  obvious  from  our  results  that  when  the  time  scale  of  basic  state  variation  is  much  less 
than  7'*,  exponential  modes  of  inertial  instability  are  largely  irrelevant. 
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Inertial  instability  can  be  expected  to  have  three  consequences  in  the  tropical  mid¬ 
dle  atmosphere:  horizontal  redistribution  of  angular  momentum  (and  potential  vorticity), 
dissipation  of  laterally  propagating  Rossby  waves,  and  vertical  mixing  of  constituents.  Be- 
catise  of  their  similarity,  local  and  global  modes  may  contribute  alike  to  mean-flow  effects, 
such  as  inertial  adjustment  (Dunkerton,  1981)  and  mixing.  (The  efiiciency  of  inertial 
adjustment  is  open  to  question.  Our  recent  numerical  results  in  the  middle  atmosphere 
demonstrate  that  an  adjustment  occurs  but  is  incomplete.  See  also  Thorpe  and  Rotunno, 
1989  for  further  discussion  in  the  context  of  frontal  CSI.)  Local  modes  will  be  more  impor¬ 
tant  for  Rossby  wave  dissipation.  According  to  Killworth  and  McIntyre  (1985),  barotropic 
processes  are  insufficient  for  time-averaged  critical  layer  absorption  as  t  —*  oo. 
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APPENDIX 


Dispersion  relation  at  complex  k 

Analysis  of  absolute  instability  in  multi-dimensional  flow  is  difiicult;  a  formal  proof  of 
absolute  equatorial  inertial  instability  using  WKB  theory  has  not  been  given  and  will  not 
be  attempted  here.  (I  doubt  whether  such  a  proof  is  possible  for  equatorial  instability 
without  approximation  to  the  governing  equations).  It  is  interesting  nevertheless  to  cal¬ 
culate  eigenfrequencies  of  (2.10)  at  complex  k.  Several  limitations  of  this  approach  should 
be  noted.  (1)  Analysis  begins  with  the  complex  dispersion  relation  defined  numerically  by 
(2.10)  with  D  =  0  rather  than  the  Laplace  transform  of  the  time-dependent  system;  this 
excludes  continuum  modes  which  may  be  important  in  reality.  (2)  The  exact  dispersion 
relation  is  determined  numerically  over  a  finite  range  of  k,  rather  than  analytically  over 
all  k.  The  cubic  approximation  of  Section  2c  (and  a  quadratic  version  of  it)  gave  similar 
results,  but  their  validity  could  not  be  guaranteed  for  all  combinations  e,  k.  (3)  The  infinite 
betarplane  could  not  be  simulated;  results  were  obtained  in  a  channel.  Channel-dependent 
modes  occur  outside  the  short-wave  cutoff  (defined  by  a  hyperbola  in  the  complex-A;  plane 
extending  outward  from  kmax  on  the  real  axis).  This  region  will  be  ignored.  (4)  Construc¬ 
tion  of  WKB  solutions  in  some  cases  may  require  matching  across  a  WKB  breakdown  point 
(Pierrehumbert,  1984);  the  full  construction  will  not  be  attempted  but  it  will  be  assumed 
that  matching  is  possible. 

To  simplify  matters,  the  shear  was  assumed  independent  of  latitude: 

rr  /  .  fl  -  V’COSx/al  ,  XX  /  X 

U^niy-Vc)  — —  =i{y-yc)'9/{x) 

and  V  s  0,  corresponding  as  closely  as  possible  to  the  theoretical  results  of  Section  2c. 
Local  stationary  instability  exists  near  ye  =  y«;  discussion  will  center  around  this  choice 
of  y«.  The  eigenproblem  (2.10)  was  solved  at  z/a  =  x,  the  location  of  maximum  shear 
(where  9/  =  1),  allowing  k  to  be  complex.  Eigenfrequencies  were  obtained  by  a  shooting 
method  for  several  combinations  7,  c.  Results  shown  m  Fig.  13  were  representative  when 
<  *«wt  (in  this  case,  7  =  1.5  •  10“*8~*,  «  =  0.1657  m“*8*).  Growth  rate  increased  in  both 
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Fig.  13:  Imriflcid  growtl'  r  '  >  as  a  fonction  of  complex  k  obtained  firom  the  eigenproblem 
(2.10)  with  P  a  C,  ^  "  'r(y  -  y,),  7  *  1.5  •  10"*8"',  €  =  0.1657  m"V,  and  y,  ss  y,. 
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directions  away  from  the  real  axis,  and  there  was  a  saddle  point  near  s  =  ka  =  (5.55,-1.50) 
-  not  far  from  the  ‘most  unstable  wavenumber’  on  the  real  axis.  Existence  and  location  of 
the  saddle  are  a  function  of  p2U’ameters  (-7,  £,  ’>£'/,  ye,  etc.)  so  there  is  nothing  special  about 
this  example.  In  the  inversion  of  Laplace  transform,  the  saddle  is  considered  ‘imavoidable’ 
when  integrating  over  kr  if  unstable  branches  originate  on  opposite  sides  of  the  real-A;  axis 
(Pierrehumbert,  1986).  This  was  apparently  the  case  based  on  numerical  results,  but  an 
analytical  proof  is  beyond  the  scope  of  this  paper. 

Not  all  saddles  are  imavoidable  in  this  problem,  particularly  those  on  the  imaginary 
axis  resulting  from  merger  of  unstable  branches  originating  in  the  upper  half-plane.  This 
sort  of  behavior  was  observed  when  ti  =  Try. 

Knowledge  of  u>(k)  is  sufficient  to  determine  new  eigenfrequencies  (j(A;)  from  a  complex 
Doppler  shift 

w(fc)  =  C>{k)  +  A:uo 

where  C)  is  given  by  the  complex  dispersion  relation  in  the  ‘control’  case  y*  =  y,  (e.g., 
Fig.  13).  Other  effects  could  be  added,  such  as  x-diffusion.  In  WKB  theory,  breakdown  of 
Galilean  invariance  is  evidently  due  to  the  imaginary  part  of  k,  as  can  be  seen  by  translating 
a  wave  packet  of  exponential  shape  to  the  left  or  right.  For  constant  shear,  variation  of 
Ve  introduces  a  constant  mean  flow  change  So  =  — iSye'9f.  The  location  of  saddle  point 
and  its  associated  eigenfrequency  change  as  a  function  of  ye;  this  is  illustrated  in  Fig.  14 
where  comparison  is  made  to  time-dependent  model  results  obtained  with  ^  =  0.25.  The 
agreement  is  rather  good.  Growth  maximized  when  the  mean  flow  was  zero  near  the  center 
of  the  unstable  region.  Due  to  noticeable  asymmetry  of  u(A;)  between  upper  and  lower 
half-planes,  results  were  not  exactly  symmetric  about  ye/y«  =  1.0. 

Growth  rates  were  slightly  smaller  in  the  time-dependent  model,  by  a  few  percent  - 
powibly  due  to  y-diffusion,  included  in  the  time-dependent  model  but  not  in  the  eigen- 
problem  (2.10).  (Horizontal  diffusivity  in  the  time-dependent  model  was  approximately 
2.1  •  10*m’s~^.  Saddle-point  values  in  Fig.  14  were  corrected  for  z-diffusion  only.)  There 
wwabo  a  non-WKB  effect  such  that  growth  rates  in  the  time-dependent  model  dimimshed 


38 


as  V'  increased,  i.e.,  as  the  size  of  the  unstable  region  wzis  reduced  in  z.  Simulations  with 
iff  >  0.25  underestimated  the  true  growth  rate  somewhat. 

Calculated  and  observed  periods  were  in  reasonable  agreement  for  ‘advected’  instabil¬ 
ities.  Note  that,  unlike  the  situation  in  Fig.  3,  Ur  was  nonzero  everywhere  except  near 
ye/y$  =  1-0;  this  was  due  to  the  smaller  value  of  'y  and  relatively  smaller  contribution  from 
the  zonally  symmetric  component  of  instability. 

Displacement  of  the  saddle  point  off  the  real  axis  implies  growth  of  the  wave  packet 
between  z/a  =  tt  and  its  center  which  is  displaced,  in  physical  space,  downstream  from  the 
region  of  most  unstable  flow.  This  expectation  was  verified  by  numerical  results,  including 
a  more  subtle  prediction  that  zonal  wavelengths  are  slightly  larger  (smaller)  for  packets 
displaced  to  the  left  (right)  of  center. 

Finally,  it  was  possible  to  vary  z/a,  holding  other  paraimeters  fixed,  and  determine 
k{u,  x)  with  w  set  equal  to  its  saddle-point  value  as  in  Fig.  10  of  Pierrehumbert  (1984).  This 
was  a  tedious  calculation  due  to  rapid  variation  of  k;  only  two  examples  were  attempted, 
one  at  CaMti  and  another  above.  In  both  cases,  it  was  striking  how  quickly  k  departed  from 
the  real  axb  in  both  directions  as  was  reduced  below  unity,  which  may  help  explain 
the  ‘locally  symmetric’  appearance  of  local  instability  having  few  cells  in  the  z-direction. 
Also,  it  may  explain  why  local  modes  can  be  foimd  in  weakly  inhomogeneous  flow,  as  only 
a  slight  reduction  of  below  unity  immediately  moves  k  into  the  opposite  half-plane  -  a 
necessary  condition  for  the  existence  of  local  modes  (Pierrehumbert,  1984). 

When  e  >>  CMuti  the  saddle-point  trajectory  is  closer  to  the  imaginary  axis,  consis¬ 
tent  with  numerical  results  demonstrating  the  preference  for  quasi-stationary,  single>cell 
instability  in  this  case. 
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